First passages in bounded domains: When is the mean first passage time meaningful? 
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We study the first passage statistics to adsorbing boundaries of a Brownian motion in bounded 
two-dimensional domains of different shapes and configurations of the adsorbing and reflecting 
boundaries. From extensive numerical analysis we obtain the probability P(u>) distribution of the 
random variable uj — ti/(ti+T2), which is a measure for how similar the first passage times ti and T2 
are of two independent realisations of a Brownian walk starting at the same location. We construct 
a chart for each domain, determining whether P(co) represents a unimodal, bell-shaped form, or a 
bimodal, M-shaped behavior. While in the former case the mean first passage time (MFPT) is a 
valid characteristic of the first passage behavior, in the latter case it is an insufficient measure for 
the process. Strikingly we find a distinct turnover between the two modes of P(co), characteristic 
for the domain shape and the respective location of absorbing and reflective boundaries. Our 
results demonstrate that large fluctuations of the first passage times may occur frequently in two- 
dimensional domains, rendering quite vague the general use of the MFPT as a robust measure of 
the actual behavior even in bounded domains, in which all moments of the first passage distribution 
exist. 

PACS numbers: 02.50.-r: 03.65.Nk: 42.25. Dd; 73.23.-b 



I. INTRODUCTION 

The concept of first passage underlies diverse stochas- 
tic processes in which it is relevant when the value of the 
random variable reaches a preset value for the first time. 
A few stray examples across disciplines include chemi- 
cal reactions [l|-[6| , the firing of a neuron 0, H| , random 
search of a mobile or immobile target |9l-[l9j. diffusional 
disease spreading [20| . DNA bubble breathing [2l| . dy- 
namics of molecular motors [22|, [23[ , the triggering of a 
stock option J24[, etc. A variety of first passage time 
phenomena and different related results have been inves- 
tigated in Rcfs. [25|, [26j]. While for continuous processes 
the first passage across a given preset value coincides with 
the first arrival to exactly this value, for Levy flights char- 
acterized by long-tailed jump length distributions with 
diverging variance both quantities become different, and 
large overshoots across a preset value occur [27j ■ 



The distribution of first passage times in unbounded 
domains is typically broad, such that not even the mean 
first passage time exists pfij . In particular, in one- 
dimensional, semi-infinite domains the first passage time 
distribution of a Markovian process is universally dom- 
inated by the i -3 / 2 scaling nailed down by the Sparre 
Andersen theorem [26j . A similar divergence of the mean 
first passage time occurs in stochastic processes charac- 
terized by scale- free distributions of waiting times [28| . 
In contrast, in many practically important situations first 
passage processes involve particles which move randomly 
in bounded domains (see, e.g., Refs. [29U321 ] ) . In this case 
the random variable of interest, the first passage time r 
to, e.g., a boundary, a target chemical group, a binding 
site on the surface of the domain or elsewhere within the 
domain, etc., has a distribution 4 r (r) of the generic, gen- 
eralized inverse Gaussian form (see e.g., the discussion in 
Ref. [l5| and references therein) 



^(t) ~ exp 
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where a and b are some constants dependent on the shape 
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of the domain, the exact starting point within the do- 
main, etc., and /x is the so-called persistence exponent 
[33j . When the linear size of the domain (say, the ra- 
dius R of a circular or a spherical domain) diverges (i.e., 
R — > 00), the parameter b also diverges such that the 
long-time asymptotic behavior of the first passage time 
distribution is of power-law form without a cutoff. In 
this case, at least some, if not all, of the moments of 
^(r) diverge. 

The first passage time distribution in Eq. (fTJ) is exact 
only in the particular case of Brownian motion on a semi- 
infinite line in presence of a bias pointing towards the 
target site, or, cquivalently, for the celebrated integratc- 
and-fire model of neuron firing by Gerstein and Mandel- 
brot Q. In general, the detailed form of ^(r) is obviously 
much more complex than given by Eq. (TTJ) , depending on 
the very shape of the domain under consideration and the 
exact boundary value problem. Typically \P(t) is given 
in terms of an infinite series. Nonetheless, on a qualita- 
tive level, the approximation (TTJ) provides a clear picture 
of the actual behavior of the first passage time distri- 
bution in bounded domains. Namely, ^(t) consists of 
three different parts: a singular decay for small values of 
r, which mirrors the fact that the first passage to some 
point starting from a distant position cannot occur in- 
stantaneously. This is followed at intermediate times by 
a generic power-law decay with exponent fj,, depending 
on the exact type of random motion. Finally, an expo- 
nential decay at long r cuts off the power-law. A crucial 
aspect is that the exponential cutoffs at both short and 
long t ensure that in bounded domains \&(t) possesses 
moments of arbitrary positive or negative order. 

Distributions of the form (TTJ) are usually considered 
narrow, as opposed to broad distributions, which do not 
possess all moments [H|-[28j], e.g., *(r) in Eq. (fT|) with b = 
00. Once all moments exist, it is often tacitly assumed 
that the first moment of this distribution, the mean first 
passage time (MFPT) 

(t) = / r*(r)dr, (2) 

is an adequate measure of the first passage behavior. The 
actual analytical calculation of the MFPT may require a 
considerable computational effort, and the calculation of 
higher moments is quite formidable and is not always 
possible, (compare, e.g., Refs. [29l - l3l"l ]). Conversely, it 
has been demonstrated in, e.g., recent Refs. [33 - l39| that 
random variables with truncated power-law distributions 
behave in several important aspects as those character- 
ized by non-truncated, broad distributions, revealing sub- 
stantial fluctuations between individual realisations and 
thus rendering the concept of a mean first passage time 
a bit unsubstantiated. To be more precise, this concerns 
not the functional form of the MFPT for a given pro- 
cess, but rather its use as a characteristic quantity for 
the process. The functional form of the MFPT is cer- 
tainly an important property, providing valuable insights 



to the scaling behavior, for instance with the system size 
or the initial distance of starting point and target. In 
contrast, the very numerical value of the MFPT can sig- 
nificantly differ from the values drawn from individual 
trajectories. Therefore, the MFPT can be substantially 
larger than the most probable value for the first passage 
time. Clearly, an understanding of how representative 
the MFPT is of the actual behavior and, concurrently, 
how important fluctuations of r between individual real- 
isations indeed are of utmost conceptual importance in 
many areas, such as, e.g., an interpretation of the first 
passage data obtained from single particle tracking. 

In this paper we analyze, via extensive Monte Carlo 
simulations the role of fluctuations between individual 
realisations of first passage times for Brownian motion 
(BM) in two-dimensional bounded domains of different 
shapes, and with different configurations of the reflective 
and adsorbing boundaries. Analogous results for three- 
dimensional systems and for systems with quenched dis- 
order will be presented elsewhere |40j . 



II. SIMULTANEITY CONCEPT OF FIRST 
PASSAGE 

To quantify the relevance of such fluctuations and the 
effective broadness of the corresponding first passage time 
distribution ^(r) we employ a novel diagnostics method 
based on the concept of simultaneity of first passage 
events, compare Fig. [TJ Instead of the original first pas- 
sage problem with quantifying the statistical outcome for 
a single Brownian walker, we simultaneously launch two 
identical, non-interacting Brownian particles at the same 
position ro (which is identical to two different realisa- 
tions of the trajectories of a single BM starting at ro). 
The corresponding outcomes are the first passage times 
Ti and T2 . We now define the random variable 

" = — (3) 

Tl + T 2 

such that u) ranges in the interval [0,1]. The unifor- 
mity index u> measures the likelihood that both walkers 
arrive to the adsorbing boundary simultaneously: when 
uj is close to 1/2, the process is uniform and the parti- 
cles behave as if they were almost performing a Prus- 
sian Gleichschritt. In contrast, values of w close to or 
1 mean highly non-uniform behavior, implying that the 
MFPT is not representative of the actual behavior, but 
is merely the first moment of an effectively broad dis- 
tribution. We note parenthetically that similar random 
variables have been used in the analysis of random prob- 
abilities induced by normalization of self-similar Levy 
processes [Il[ , of the fractal characterization of Paretian 
Poisson processes [13] , and of the so-called Matchmaking 
paradox [13, El- 

Within a given bounded domain we evaluate the distri- 
bution P(lj) measuring the uniformity of the first passage 
dynamics with respect to some fixed starting point Tq. 
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FIG. 1. (Color online) Trajectories of two Brownian walkers 
starting at the same initial position (ro,#o) inside a bounded 
pie-wedge domain with opening angle © as well as absorb- 
ing radial boundaries (dashed lines) and reflecting boundary 
(solid line) at r — 1. The values of the first passage times to 
the adsorbing boundaries are used to construct the random 
variable u. 

This is repeated for a large number of nodes ro within 
the domain, thus producing a uniformity chart of first 
passage. Remarkably, we find that the very shape of this 
distribution depends delicately on the domain shape, the 
actual settings of adsorbing and reflecting boundaries, 
and on the starting location ro. In some starting ar- 
eas P(u>) has a characteristic unimodal, bell-shaped form 
with a maximum at w = 1/2, signaling that most pairs 
of BMs will arrive to the adsorbing boundary simultane- 
ously. This means, in turn, that in this case the parental 
first passage time distribution ^(r) can be considered as 
sufficiently narrow such that the MFPT can be consid- 
ered as a plausible measure of individual first passage 
events, providing a rather accurate estimate for the typ- 
ical value of the first passage time. Conversely, we find 
that for other starting areas P(ui) exhibits a completely 
different behavior and has a characteristic bimodal, Mi- 
shaped form with a local minimum at a; = 1/2 and two 
maxima close to uj = and to = 1. In that case simul- 
taneous arrival of two initially synchronized walkers is 
unlikely, i.e., any two trajectories will most likely possess 
distinctly different first passage times. The parental first 
passage time distribution \&(r) is consequently broad and 
sample-to-sample fluctuations matter: the MFPT can- 
not be considered as an adequate measure of the actual 
behavior. Given that, by definition, the averages (ti) 
and (T2) are identical and, moreover, the moments of t\ 
and T2 of arbitrary order coincide, one can think of to 
(and, hence, of the distribution P(w)) as a measure of 
the symmetry breaking between different realisations of 
the process. Note also that situations in which the mean 
value of some pertinent parameter is dominated by the 
tails of the distribution, and this mean thus has a very 
different value compared to the most probable value (and 
may even show a completely different dependence on the 
system parameters) is most often encountered in disor- 



dered systems [H, • Here we observe such a behavior 
in absence of any disorder. 

Scanning then over the possible starting points within 
each bounded domain, we obtain a corresponding phase- 
chart for P(w), distinguishing regions in which P(uj) 
has M-shaped or bell-shaped behavior. The demarca- 
tion zone between these two phases, depicted by beige 
color in the following, represents a plateau-like, almost 
uniform behavior of P(oj) with zero second derivative at 

LO = 1/2. 

We proceed by giving a general definition of the first 
passage time distribution 4 r (r) and its corresponding 
MFPT in Section IIIIl and also establish a relation be- 
tween ^(t) and the uniformity distribution P(cj). In 
Section IIVI we study in detail the problem of Brownian 
motion in a pie- wedge shaped domain with absorbing and 
reflecting boundaries. In Sections [V] and I VII we discuss 
the forms of P(w), as a function of the location of the 
starting point, for circular domains with small aperture 
on the boundary, a two-dimensional version of the so- 
called Narrow Escape Time problem [48}, and for trian- 
gular domains with adsorbing boundaries, respectively. 
Our results are summarized in Section IVIII 



III. FIRST PASSAGE DISTRIBUTION, MEAN 
FIRST PASSAGE TIME, AND THE 
UNIFORMITY DISTRIBUTION P(uj). 

Consider a BM inside a general two-dimensional do- 
main <S, whose boundary dS = dS a U dS r comprises re- 
flecting, dS r , and absorbing, dS a , parts. At time t = 0, 
the BM initiates at ro £ 5 and evolves within the do- 
main until the trajectory hits dS a for the first time at 
some random instant r. Furthermore let P(r, £|ro) de- 
note the conditional probability distribution for finding 
the Brownian walker at position r at time t, provided 
the initial condition was at ro at t = 0. The distribution 
P(r, i|r ) is the solution of the diffusion equation 

^P(r,£|r ) = ^V2p(r,i|r ) (4) 

on iS, where V 2 is the two-dimensional Laplacian equiva- 
lent to d 2 /dx 2 +d 2 /dy 2 in Cartesian coordinates. Eq. Q 
is subject to the initial condition as well as the bound- 
ary conditions at dS. Here D is the diffusion coefficient. 
The solution of this boundary value problem is, in the 
best case, cumbersome, and explicit solutions may be ob- 
tained for only few simple geometries, compare Ref. [45| . 

If a finite part of the boundary is absorbing, i.e., dS a 
is not empty, then the distribution P(r, i|r ) is no longer 
normalized. The survival probability 5^ TQ (t) that the 
walker has not reached dS a up to time t, is defined by 

^r (t) = f P(r,t\r )dv. (5) 
Js 

5^ ro [t) is a monotonically decreasing function of time, 
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eventually reaching zero value, lim t _ ) . 00 o$" ro (i) = 0. The 
desired distribution of first passage times to the adsorb- 
ing boundary becomes 



*r (r) = - 



d^ro(r) 
dr ' 



(6) 



The MFPT associated with the distribution *(t) is de- 
fined as the first moment 



(r)(r ) 



T^ ro (r)dT 



y ro (r)dr. (7) 



We note parenthetically that in most of the existing lit- 
erature, apart of recent Refs.pl, l29l. l3(i|. the dependence 
of the MFPT on the starting position of the walker is ei- 
ther simply neglected, or it is assumed that the starting 
point is randomly distributed within the domain S. As 
we proceed to show, the rvdependence of the first pas- 
sage time distribution is a crucial aspect which cannot 
be neglected. 

We now turn to the uniformity distribution P(lo) of 
the random variable u, Eq. Let 



$(A) = / P(u) exp (-Aw) dui, 
Jo 



(8) 



with A > 0, denote the moment generating function of uj. 
Since n and r 2 are independent, identically distributed 
random variables, expression © can formally be repre- 
sented as 



where K^^iS) is the modified Bessel function of the sec- 
ond type. It was realized [HI, HH that the form of the 
distribution P(oj) in Eq. (fT2| is distinctly sensitive to the 
value of the persistence exponent p, which characterizes 
the scaling behavior of the first passage time distribu- 
tion \&(t) at intermediate times. Thus, for p > 1, P(uj) 
is always a unimodal, bell-shaped function with a max- 
imum at uj = 1/2. For p = 1, P(uj) is almost uniform, 
P(uj) ~ 1, apart from narrow regions at the corners uj = 
and uj = 1, for b/a ^> 1. Curiously, for p <1, which cor- 
responds to the most common case, there exists a critical 
value p c of the ratio p = b/a such that for p > p c the dis- 
tribution P(oj) has a characteristic M-shapcd form with 
two maxima close to and 1, while at w = 1/2 we find 
a local minimum. Such a transition from a unimodal, 
bell-shaped to bimodal, M-shaped form mirrors a signifi- 
cant manifestation of samplc-to-samplc fluctuations that 
has been indeed observed in exact calculations of P(uj) 
for Brownian search processes for an immobile target in 
d-dimensional spherical geometries [l5[ . 

In what follows we further explore this intriguing be- 
havior of the first passage time distribution via extensive 
Monte Carlo simulations focusing on the effects of the 
domain shape, the type of the boundary conditions, and 
the initial position of the walker. 



IV. UNIFORMITY DISTRIBUTION P(u>) IN A 
PIE- WEDGE DOMAIN 



$(A)= / / *(r 1 )*(r 2 )exp(-A^^ 

Jo Jo \ n + t-2 



d Tl dT 2 . (9) 



Integrating over dr\ we change the integration variable, 
Ti — > uj, so that Eq. ((SJ) is rewritten in the form 



$(A) = J exp (-Aw) — 



du> 



t 2 *(t 2 )* 



l-OJ 



dro 



(10) 



From comparison with Eq. ([5]), we readily read off the 
desired distribution function 



P(oj) 



1 



dr. (11) 



Therefore, P(oj) is known for given ty(t). 

To get an idea of the typical behavior of the uniformity 
distribution P(oj), we use the generic form ([TJ for the first 
passage time distribution. From Eq. (|11[) we find from 
integration that 



P(oj) 



2KI 2Vf w(l - uj) 



(12) 



Consider now the case of a bounded domain of pie- 
wedge shape with unit radius, R = 1 and opening angle 
G. The absorbing boundaries correspond to the radial 
edges, while the outer circular edge is reflective, compare 
Fig. [1] Clearly, for a BM inside such a pie- wedge domain, 
all moments of the first passage time distribution exist. 

Before we proceed to investigate this case we first turn 
to the case when the wedge radius is infinite, R — > oo. 
Then the distribution function P(r, i|ro) is known ex- 
actly, (see, e.g., Refs. [H, H^|) and is represented by an 
infinite series whose leading term for t — > oo is given, up 
to a normalization constant, by 



P(r,t\r ) ~ 



Trsin {ttOq/Q) tf+ftytpt 



ADQt 
( PoP\ 
\2Dt) 



(13) 



where I v (z) is the modified Bessel function of the first 
kind and r = (p, 6) is conveniently represented in polar 
coordinates. This solution is obtained for the sharp ini- 
tial condition P(r, 0|ro) = 7rsin(7r#o/0)^(r — ro)/26po- 
From Eq. (fTB"]) one finds the asymptotic behavior of the 
survival probability, 



TT/2B 



t 



-tt/20 



(14) 
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such that the first passage time distribution becomes 

(15) 

Note that this distribution is of the generic form (pj, 
where b = oo due to the infinite domain size. The non- 
universal persistence exponent is given by 



Therefore, the MFPT diverges when O > tt/2 and is 
finite for 6 < tt/2. According to the qualitative analysis 
from Section IIII1 P(uj) will have a bimodal form in the 
former case and a unimodal one in the latter. 

We now turn our attention to finite-sized pie-wedges, 
for which the MFPT and all higher moments of Eq. © 
are finite. In principle, an exact solution for the first 
passage time distribution in this case can be obtained 
by solution of the corresponding mixed boundary value 
problem, but the result will be too cumbersome for our 
purposes. Instead, we resort to numerical simulations. 
We now show that for finite pie-wedges the actual be- 
havior is in fact richer than in the case of an infinite 
wedge. 

We performed Monte Carlo simulations of a random 
walk inside a pie- wedge of unit radius and opening angle 
0. The boundary conditions along the radii are absorb- 
ing and reflecting along the circular edge, see Fig.Q] The 
random walk is simulated in terms of a standard Pearson 
walk on a plane (compare Ref. |47j]), which consists of a 
sequence of steps of fixed length A = 0.001 and uniform 
waiting time v = 1/A. After each step the walker turns 
by a random angle with uniform distribution. At time 
t = the walker is released at (po, Oo), and its trajectory 
is recorded until it hits a point on the absorbing bound- 
ary for the first time. Generating N (we used N = 10 5 ) 
such trajectories, we obtain a set of first passage times 
{Ti}, from which we construct the first passage time dis- 
tribution. Since all n arc independent, identically dis- 
tributed random variables, the uniformity distribution 
P(ll>) is then readily obtained via Eq. (|3J) from distinct 
pairs Ti and T2 chosen at random from the set {t^}. 

Fig. |5fa) shows the first passage time distributions cor- 
responding to a fixed po = 0.76 and three different start- 
ing angles 9q for a pie- wedge with opening angle = tt/2. 
One notices that, for small and large values of r, all three 
distributions &(t) significantly deviate from the inter- 
mediate power-law behavior, which is due to exponential 
tempering. On the other hand, at intermediate times 
the distributions exhibit a slower, power-law like decay 
within a range that depends significantly on 6q. The nar- 
rowest distribution (light blue) is obtained for a starting 
position (0.76, 0), which is exactly on the symmetry axis 
of the wedge. Increasing the angle 9o away from the 
symmetry axis results in a broadening of &(t), and the 
intermediate algebraic decay is more pronounced. 

In panel (b) of Fig. [5] we plot the corresponding uni- 




FIG. 2. (Color online) (a) First passage time distribution 
^( P0 ,e )( T ) f° r a Brownian walk in a pie-wedge domain with 
opening angle = tt/2. Different colors and symbols corre- 
spond to different starting points: (po — 0.76, 9q — —0.38 
(dark blue/triangles), (0.76,-0.23 0) (beige/squares), and 
(0.76,0) (light blue/circles). The dashed straight line in- 
dicates the intermediate power-law decay *&(t) ~ V 1 " 2 , 
Eq. (|15|) . (b). The corresponding distribution P(lj) with 
the same color and symbol coding. In grey scale, dark blue 
corresponds to the darkest shade of grey, light blue to dark 
grey and beige to light grey. 



formity distributions P(uS). For the narrowest first pas- 
sage time distribution ^(r) (light blue symbols), P(uj) 
is bell-shaped with 1/2 representing the most probable 
value, such that sample-to-sample fluctuations of r are 
less significant. In this case, apparently, the MFPT is a 
meaningful, reliable measure of the first passage behav- 
ior, and any two walkers starting from the same position 
on the symmetry axis of the wedge will most likely be 
absorbed at the same instant of time. Strikingly, we find 
that this is no longer valid for the two other starting po- 
sitions off the symmetry axis: for ro = (0.76, — tt/3) the 
uniformity distribution P(ui) is almost uniform, except 
for narrow regions in the vicinity of the edges, meaning 
that any relation between the first passage times of two 
walkers is equally probable. Finally, for starting posi- 
tion (0.76, — 7r/7), which is the one closest to the absorb- 
ing boundary P(ui) has a characteristic M-shaped form 
with maxima close to and 1, and a local minimum at 
uj = 1/2. This signifies that in a such a case the sym- 
metry between any two walkers is most distinctly broken, 
and they will be absorbed at very different times. Clearly, 
in this case the MFPT is not representative of the actual 
behavior. 

To quantify the shape of the distribution of uj we per- 
form a fit of the numerically obtained P(u>) to a quadratic 
polynomial of u> in the domain 0.05 < cj < 0.95. From 
this fit we obtain the coefficient x of t ne quadratic term. 
The sign of x thus determines the shape of P(u>): x < 
corresponds to the unimodal, bell-shaped distribution, 
X > signifies that the distribution is bimodal, M- 
shaped, and a zero value of x — means that P(uS) is 
uniform. Such a procedure, of course, has some ambigu- 
ities, especially when we deal with the demarcation line 
X = between regions in which P(u>) has unimodal and 
bimodal forms, as it is not always clear how many digits 
arc to be taken into account. This results in a certain 
broadening of the demarcation line. However, we have 
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FIG. 3. Dependence of the parameter \ on the starting angle 
6*o for fixed po = 0.76 in a pie-wedge domain with opening 
angle = n/2. The three insets show the shape of the uni- 
formity distribution P{fJ) for the 9o values indicated by the 
arrows. 




FIG. 4. (Color online) Phase-chart for the shape of the unifor- 
mity distribution P(lj) in three different pie-wedge domains: 
(a) 9 = 3tt/4, (b) 6 = tt/2, (c) 9 = tt/3. The color of the 
symbols is light blue if \ < — \*, dark blue if x > X*i an d 
beige if \x\ < X+> where we chose x* = 0.25. For each ini- 
tial location, -P(oj) was computed from a sample of N — 10° 
random trajectories. 



checked in several cases that this procedure produces re- 
liable results. In Fig. [3] we show the evolution of \ versus 
the starting angle 9q for fixed po = 0.76 in a pie-wcdgc 
domain with the opening angle O = tt/2. We observe a 
continuous, periodic variation of P(uj) with the starting 
position, changing from an M-shaped to a bell-shaped 
form. The insets of Fig. [3] show the schematic distribu- 
tion P(u>) for some specific values of 8q. The absolute 
value of x indicates how far the distribution P(lu) devi- 
ates from a locally uniform distribution. 

Finally, we used this approach to create the phase- 
chart for the shape of the uniformity distribution P(co) 
with respect to the starting position of the walker within 
the pie- wedge domain. In Fig. [4] we present a systematic 
scan of the domain for three pie-wedges with different 
opening angles. One observes that in all three cases, 
there exists a region in which P(ui) is bell-shaped (light 
blue symbols) and a region with M-shaped P(w) (dark 
blue symbols), separated by a small region with nearly 
uniform distribution (beige symbols). 

Therefore, as we have already remarked, the actual 
behavior in a finite pie-wedge appears to be much richer 
than in an infinite wedge. Consider an experiment in 
which one aims to find an estimate of the MFPT by 



tracking the evolution of a few single particle trajecto- 
ries starting at the same position inside the light blue 
region. The outcome of such an experiment will be a 
good estimate of the MFPT, with reliably small error. 
This will be the case since in the light blue region P(io) 
is bell-shaped, which means that the probability that the 
two trajectories arrive at the same time is maximal. In 
contrast, if two single particle trajectories start anywhere 
inside the dark blue region, then it is most likely that 
these trajectories will arrive to the adsorbing boundary 
at very different times, yielding a poor and unreliable es- 
timate for the MFPT. The sample-to-sample fluctuations 
in this case are very important and, as a consequence, the 
MFPT is not an adequate measure of the actual behav- 
ior. Qualitatively, the sample-to-sample fluctuations of 
the MFPT increase as the trajectories start closer to the 
absorbing boundaries. However, this is not always true, 
as can be observed in Fig. 2] (c) for the pie-wedge with 
= 7r/3 for which the light blue region extends toward 
the vertex of the wedge. 



V. CIRCULAR DOMAIN WITH APERTURE 

We now turn our attention to the first passage time 
problem of a Brownian particle in a circular domain of 
unit radius and the following boundary conditions: the 
segment with \9\ < 0/2 is absorbing while the remaining 
part of the outer circle is reflective. The aperture of the 
circular domain is thus of angle 0. One often encounters 
a three-dimensional version of this problem in cellular 
biochemistry, when one is interested in the time needed 
for a particle (a ligand, etc.), diffusing within a bounded 
domain (for instance, a microvesicle) to reach a small es- 
cape window or a binding site, which is an aperture in 
an otherwise reflecting boundary. This is the so-called 
Narrow Escape Time problem, which attracted consid- 
erable attention within the last two decades (see, e.g., 
Refs.[3ll. |48| and references therein). 

We analyze the shape of the uniformity distribution 
P{uj) as a function of the starting point of a Brownian 
walker. As in the previous section, we generate N = 10 
random walks commencing from the same starting posi- 
tion (po,9o) inside the unit circle and determine the set 
{Tj} of first passage times to the location of the aper- 
ture. From these data we obtain P{u>) and compute 
the parameter x- I n Fig. [S]wc show the phase-chart for 
the shape of P(ui) for four different sizes of the aper- 
ture: = 7r/18, 7r/2, 7t, and 2ir. Each symbol in the 
charts is light blue, beige, or dark blue, depending on 
whether the corresponding starting position leads to a 
bell-shaped, uniform, or M-shapcd distribution P(uj), re- 
spectively. Note that the case shown in Fig. [5] (d) reduces 
to a one-dimensional problem (see, e.g., Rcf. |15|). 

Similarly to our findings for the pie-wedge domain, 
we observe that the MFPT is not always a represen- 
tative measure for the two-dimensional Narrow Escape 
Time problem. Interestingly, the sub-domain in which 
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FIG. 5. (Color online) Phase-chart for the shape of the unifor- 
mity distribution P(oj) for a Brownian walker in the unit circle 
with reflective BCs (solid lines) and an aperture of size (ab- 
sorbing BC - dashed lines), (a) 9 = tt/18, (b) 6 = tt/2, (c) 
O = 7r, and (d) = 2n. Starting locations are colored light 
blue if x < ~~ X*j dark blue if \ > X* an d beige if \x\ < X*i 
where \* = 1- in the lower panel, the relative error of the 
FPT e is shown as function of the initial radial position r 
along the horizontal diameter of the unit circle with = 7r/2, 
panel (b): r = is associated to the center of the circle and 
the point r = 1 lies over the absorbing boundary. 



the MFPT is the least probable outcome (dark blue cod- 
ing) is practically the same for small holes of 9 < w/2. It 
is worthwhile noting that while in this region P(ui) is al- 
ways bimodal, the height of its maxima increases (and its 
value P{oj = 1/2) representative of the MFPT decreases) 
depending on the distance from the opening. 

In addition, from the set of first passage times {rj} we 
directly computed the MFPT and the variance vai(r). 
Both statistical indicators grow with the distance from 
the absorbing boundary. A more sensitive measure is the 
relative error e, defined as the ratio 



In the lower panel of Fig. [5] we show the dependence of 
the relative error on the starting position ro for trajec- 
tories initiating along the symmetry axis of the domain, 
namely with respect to r for fixed 9 = 0. In agreement 
with the qualitative results of the phase-chart, e < 1 only 
when P(oj) is bell-shaped, and the MFPT is the most 



FIG. 6. (Color online) Phase-chart of the shape of the unifor- 
mity distribution P(uj) for symmetric triangles with central 
angle n/2 (a) and 2tt/3 (c), and for the asymmetric triangle 
with angles 2tt/3, 7r/4, 7r/12 (b). The color of the symbols 
is light blue if \ < — X*> dark blue if x > X*i an d beige if 
|x| < X*, with X * = 0.25. 

probable outcome of a single-particle trajectory, namely, 
for trajectories starting far enough from the absorbing 
boundary. Clearly, the closer the starting position is to 
the absorbing boundary the larger the relative error be- 
comes. Very near the absorbing boundary the standard 
deviation of the first passage time becomes much larger 
than its mean. We note that this result is generic irre- 
spectively of the aperture size. 

However s is just a number and it is not clear how 
to interpret it. For instance, e = 1 or e = 2, are these 
values too small or large enough to allow us to say that 
trajectory-to-trajectory fluctuations are significant? On 
the other hand, the distribution of the simultaneity index 
which we discuss here gives a lucid answer on this ques- 
tion as manifested by the change of modality of P(oj). 

VI. TRIANGULAR DOMAIN WITH 
ABSORBING BOUNDARIES 

Finally, as a complementary example we consider a 
domain whose boundaries are completely absorbing. We 
consider the triangular domains shown in Fig. [5J two 
symmetric triangles with central angle 7r/2, panel (a), 
and 27r/3, in panel (c), and an asymmetric triangle with 
angles 27r/3, 7r/4, 7r/12 shown in panel (b). 

In Fig.[5]we show the phase-chart of the shape of P(u>), 
the results obtained arc qualitatively the same as for the 
previous two examples with mixed boundary conditions. 
P(u>) is bell-shaped and the MFPT the most probable 
outcome only when the trajectory starts far enough from 
the boundary. Also, similarly to what we observed for the 
pie- wedge domain, the domain in which P(uj) is unimodal 
extends toward the absorbing boundary if the vertex an- 
gle is less than 7r/2. 

VII. CONCLUSIONS 

We explored the problem of first passage of a Brow- 
nian particle to the absorbing boundary of finite, two- 
dimensional domains. From our study of the character- 
istic shapes of the associated distribution of the unifor- 
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mity index uj we demonstrated that the MFPT represents 
the most probable outcome (and thus is quite meaning- 
ful) only if the trajectories start in a certain subregion 
of the total domain. For starting points in the comple- 
mentary region the MFPT becomes the least probable 
outcome, indicating very large sample-to-sample fluctu- 
ations. These observations are generically important for 
single trajectory analysis of first passage time processes. 

We showed that the associated separation into bell- 
shaped and M-shaped forms of the uniformity distribu- 
tion P(ui) is a robust property of Brownian motion by 
studying the problem in different symmetric and asym- 
metric domains with mixed or fully absorbing bound- 
aries. We found that in general, sample-to-samplc fluc- 
tuations of the first passage time increase when the tra- 
jectories start close to the target boundary, leading to 
the unexpected conclusion that in such situations the 
MFPT yields insufficient information, particularly, if the 
absorption time is extracted from the outcome of very 
few single-particle trajectories. 

Next, it is worthwhile mentioning that in many inter- 
esting situations the starting position of the trajectories 
are randomly distributed inside the finite domain. From 
such analysis the so-called global MFPT is usually de- 
rived, see, e.g., Ref. [3l|. Here we found that averaging 
the associated uniformity distribution P(lu) over the do- 
main, 

P av H = / P ri »dr , (18) 
Js 

attains a uniform shape, except near ui = and lu = 1. 
This appears to be a general property of P av (w) associ- 
ated with the probability conservation, and leads to the 
unexpected conclusion that the global MFPT has little 



meaning in such situations. 

As a final remark, we emphasize that the approach 
outlined here is not limited to first passage phenomena 
only, but can be quite generally applied to probe the 
significance of sample-to-sample fluctuations of arbitrary 
random variables having distributions for which all mo- 
ments exist. Such distributions, as shown in our work, 
may appear w-broad, in the sense that the correspond- 
ing uniformity distribution P(ui) is bimodal, or, alterna- 
tively, w-narrow with unimodal P(u)). We recall that the 
variable ui has a very lucid physical meaning and its dis- 
tribution can be determined if the parental distribution 
of the random variable is known. Indeed such sort of het- 
erogeneity analysis has very recently started to be used to 
quantify sample-to-sample fluctuations in mathematical 
finances (3(| [37| , chaotic systems [H| , analysis of distri- 
butions of the diffusion coefficient of proteins diffusing 
along DNAs [H and FPT phenomena [Hj|. 
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